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ABSTRACT 

We present new high-speed, multi-observatory, multi-instrument photometry of 
the eclipsing polar UZ For in order to measure precise mid-eclipse times with the aim 
of detecting any orbital period variations. When combined with published eclipse times 
and archival data spanning ^27 years, we detect departures from a linear and quadratic 
trend of ~60 s. The departures are strongly suggestive of two cyclic variations of 16(3) 
and 5.25(25) years. The two favoured mechanisms to drive the periodicities are either 
two giant extrasolar planets as companions to the binary (with minimum masses of 
6.3(1.5)Mj up and 7.7(1. 2)Mj up ) or a magnetic cycle mechanism (e.g. Applegate's 
mechanism) of the secondary star. Applegate's mechanism would require the entire 
radiant energy output of the secondary and would therefore seem to be the least likely 
of the two, barring any further refinements in the effect of magnetic fields (e.g. those 
of Lanza et al.). The two planet model can provide realistic solutions but it does 
not quite capture all of the eclipse times measurements. A highly eccentric orbit for 
the outer planet would fit the data nicely, but we find that such a solution would be 
unstable. It is also possible that the periodicities are driven by some combination of 
both mechanisms. Further observations of this system are encouraged. 

Key words: accretion, accretion discs - methods: analytical - binaries: close - novae, 
cataclysmic variables - X-rays: stars, planetary systems. 



1 INTRODUCTION 

Approximately 20% of the known cataclysmic variables 
(CVs, see the catalogue of Ritter & Kolb 2003) are polars, 
where the primary white dwarf has a sufficiently strong mag- 
netic field to lock the system into synchronous rotation with 
the red dwarf secondary and to prevent completely the for- 
mation of an accretion disc. 

The material from the secondary overflowing the Roche 
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lobe initially falls towards the white dwarf following a ballis- 
tic trajectory until, at some distance from the white dwarf, 
the magnetic pressure overwhelms the ram pressure of the 
ballistic stream. From this point on the accretion flow is con- 
fined to follow the magnetic field lines of the white dwarf. 
The now supersonic accreting material suddenly becomes 
sub-sonic at a shock region, which forms at some height 
above the white dwarf surface. The shock-heated material 
reaches temperatures of ~ 10 — 50 keV and is therefore 
ionised. The hot plasma cools by X-ray cooling, in the form 
of bremsstrahlung radiation. With sufficiently strong mag- 
netic fields we find also cyclotron cooling, in the form of op- 
tical/infrared cyclotron radiation (e.g. ST LMi: Imamura, 
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Steiman-Cameron & Wolff, 2000, and Campbell, Harrison, 
Mason, Howell & Schwope, 2008). See e.g. Warner (1995) for 
a review of CVs and Cropper (1990) and Patterson (1994) 
for reviews of magnetic CVs. 

UZ For is one of 15 known eclipsing polars and was dis- 
covered with EXOSAT (EXO 033319-2554.2) as a serendipi- 
tous X-ray source (Giommi et al. 1987; Osborne et al. 1988). 
Extensive followup observations at multi- wavelengths estab- 
lished a ~126.5-min orbital period, one or two accretion 
spots depending on accretion state, with magnetic fields of 
~53 MG and ~48 MG and a white dwarf mass of ~0.7M Q 
(Bcuermann, Thomas & Schwope 1988; Berriman & Smith 
1988; Ferrario et al. 1989; Bailey & Cropper 1991). The 
eclipses of the accretion spots in UZ For are particularly 
rapid at 1-3 s and can only be resolved with high speed pho- 
tometry. Bailey & Cropper (1991) were able to resolve the 
eclipse of the white dwarf photosphere during a low state 
and Perryman et al. (2001) were able to resolve two accre- 
tion spots during a higher accretion state. These distinct 
rapid photometric transitions are ideal for making accu- 
rate timing measurements and therefore searching for any 
long term period variations in UZ For. Some of the above- 
mentioned authors have combined their observations with 
previous eclipse measurements in order to obtain accurate 
eclipse ephemerides. In general, significant residuals were 
seen in the O-Cs (Observed - Calculated) of the orbital pe- 
riod, but no overall trend had been detected (e.g. Perryman 
et al. 2001). More recently, Dai et al. (2010) claim the ex- 
istence of a third body orbiting UZ For in order to explain 
the O-C. However, their singular new eclipse measurement 
and subsequent derived orbital parameters are grossly in- 
compatible with all of our new observations spanning ~10 
years. 

Nevertheless, recent results of long term studies of some 
CV related objects are beginning to show trends. Parsons 
et al. (2010) presented high-speed ULTRACAM photome- 
try of 8 post-common-envelope-binaries. They detect signif- 
icant departures from linearity in some of these systems and 
suggest magnetic braking or a third body as possible mech- 
anisms to drive the O-Cs. High precision eclipse measure- 
ments of the post-common envelope binary NN Ser (Beuer- 
mann et al. 2010a) shows strong evidence for two additional 
bodies superposed on the binary's linear ephemeris. Signif- 
icant and complicated departures from a linear ephemeris 
have also been seen in the eclipsing polar HU Aqr (Schwarz 
et al. 2009). They find that neither a sinusoidal nor a 
quadratic ephemeris are sufficient to describe their O-C de- 
partures, thus more eclipse observations over the next few 
years will be needed in order to refine the ephemerides. Qian 
et al. (2010) discovered that the O-C curve of the eclipsing 
polar DP Leo shows a cyclic variation with a period of 23.8 
years. They claim that this is as a result of a giant extrasolar 
planet orbiting DP Leo, recently refined by Beuermann et 
al (2010b). 

Here we present new high-speed HIPPO, BVIT, 
SALTICAM and UCTPOL photometry of UZ For, spanning 
10 years, and use these observations to determine accurate 
mid-eclipse times of the main accretion spot in UZ For. We 
combine these with previous mid-eclipse times, that we ei- 
ther measure from archival data or extract from the litera- 
ture, and analyse for any period variations in UZ For. Note 



that adding these newer data to timing from the literature 
gives a baseline of 27 years. 



2 OBSERVATIONS 

All of the eclipse times extracted from the literature were 
published as Heliocentric Julian Dates (HJD). We have as- 
sumed that the Coordinated Universal Time (UTC) system 
was used in all cases as this was not explicitly stated in any 
of the publications. We re-corrected all times for the light 
travel time to the barycenter of the solar system, converted 
to the barycentric dynamical time system (TDB) and the 
times are listed (table 1) as Barycentric Julian Date (BJD; 
see Eastman, Siverd & Gaudi 2010 for achieving accurate 
absolute times and time standards). By doing this we have 
removed any timing systematics, particularly due to the un- 
predictable accumulation of leap seconds with UTC, and 
effects due to the influence of primarily Jupiter and Saturn 
when heliocentric corrections only are applied. We either 
calculate or re-calculate appropriate errors depending on the 
S/N and time resolution at the time of the spot ingress and 
egress. Table 1 also lists the eclipse width of the accretion 
spot and the observatory /instrument used. 

All of our new observations were also converted to 
BJD. We note that our new ground based instruments were 
synchronised to GPS to better than a milli-second. Given 
the high-speed nature of these instruments, their timing 
accuracies have been verified through simultaneous multi- 
instrument observations. The remaining space observatories 
have documented reports on the performance of their on- 
board clocks. 

2.1 Eclipse times from the literature 

The earliest UZ For eclipse measurements were made us- 
ing EXOSAT and published by Osborne et al. (1998). The 
data are of poor time resolution but are at a sufficiently 
early epoch to provide constraints for model fitting. Beuer- 
mann et al. (1988) and Ferrario et al. (1989) observed mul- 
tiple eclipses spectrophotometrically. These are also of very 
low time resolution, however the combination of multiple 
eclipses provides usable data. Allen et al. (1989) presented 
the first high speed photometry that could resolve the ac- 
cretion spot. A typographical error in the eclipse time was 
corrected by Ramsay (1994). This was soon followed with 
more high quality optical low-state photometry by Bailey & 
Cropper (1991) and high-state photometry by Imamura & 
Steiman-Cameron (1998) and EUVE light curves by War- 
ren et al. (1995). An additional high quality ST J eclipse was 
also published by de Bruijnc (2002) and three more, with 
the same instrument by Perryman et al. (2001). From the 
latter data we were able to obtain the raw observations and 
re-measure and confirm the eclipse times. 

2.2 ROSAT (1991) 

Observations were retrieved from the HEASARC archive 
and events were extracted using an aperture centered on 
the source and also a background region. The resulting light 
curves were subtracted after appropriate scaling from the 
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differing areas. The observations spanned ~1 day and mul- 
tiple eclipses were covered. The reduced data were binned 
and folded on the orbital period in order to produce a single 
eclipse light curve from which the eclipse was measured. 

2.3 HST FOS (1992) 

UZ For was observed by HST with FOS on 11th June 1992 
in two consecutive, RAPID mode observations consisting of 
925 spectra each with 1.64s exposure times. This data set 
was originally published in Stockman & Schmidt (1996), but 
their mid-eclipse times were not quoted. Therefore the HST 
data products from this observation were downloaded from 
the HST archive at the MAST. The flux- and wavelength- 
calibrated individual spectra were extracted and the flux 
summed-up between 1255 and 1518A(far UV) to create a 
lightcurve. The far UV part of the spectrum was chosen as 
this seems to have the least contribution from the accretion 
stream. The end-times for each spectra were obtained from 
the observation header keywords and the group-delay-time 
subtracted in order to obtain times of start of each expo- 
sure. Two consecutive eclipses were observed and folded and 
binned into a single eclipse from which measurements were 
made. 



2.4 EUVE 1993 and 1995 

UZ For was observed with EUVE on the 18th Novem- 
ber 1993 and on the 15th Jan 1995 for 102ks and 
76ks respectively. These data were retrieved from the 
STScI archive and reduced following the recipe from 
http://archive.stsci.edu/euve/. Lightcurves were produced 
using the xray.xtiming package in IRAF with a Is time res- 
olution. For each observation, multiple eclipses were covered 
which were folded and binned into two single eclipses from 
which measurements were made. 



2.5 UCTPOL 2002 and 2005 

One unfiltered and two BG39 filtered eclipses were obtained 
in 2002 with the University of Cape Town photo-polarimeter 
(UCTPOL) at 10 second time resolution. Two eclipse times 
were extracted: one from the unfiltered and the second from 
the folded BG39 filtered eclipses. Three unfiltered eclipses 
were obtained in 2005 at 10 and 1 second time resolution. 
Two eclipse times were extracted. On both occasions, simul- 
taneous linear and circular polarimetric observations were 
also made and reduced as in Cropper (1985). 

2.6 XMMOM 2002 

Observations were made in fast-mode using the XMM- 
Newton Optical Monitor (Mason et al 2001) between Aug 7 
and Aug 8 2002; two orbits after UCTPOL observations. The 
UVW1 filter was used (center wavelength 2910 A, FWHM 
500 A) and the data were reduced using omf chain running 
under SAS v9.0. Although three consecutive eclipses were 
observed, two had incomplete coverage. Nevertheless, upon 
folding and binning, a high signal-to-noise eclipse profile was 
obtained and eclipse measurements taken. 



2.7 SWIFT 2005 

Observations were made in event mode using the SWIFT 
UV Optical Telescope (Roming et al 2005) between Feb 2 
and Feb 6 2005. There were a number pointings, some lasting 
a few 100 sec and others a few 1000 sec. The U filter (cen- 
ter wavelength 3450A, FWHM 875A) and the V filter were 
used. Light curves were extracted using apertures centered 
on the source (radius 3 ) and also a source free background 
region with much larger aperture radius. The light curve was 
generated by suitably scaling the size of the apertures. Two 
eclipses were observed in full, one of which was simultaneous 
with the UCTPOL observations. 



2.8 SALTICAM 2007 

UZ For was observed with SALTICAM (O'Donoghue et al. 
2006) on SALT on 12 November 2007. SALTICAM was in 
slot-mode configuration, allowing a time resolution of Is 
with no deadtime. The data were reduced using the SALT 
slottools data reduction package (Crawford et al. 2010). One 
eclipse of high time resolution and signal to noise was ob- 
served from which measurements were made. 



2.9 BVIT 2009 

BVIT (Berkeley Visible Imaging Tube: Siegmund et al. 
2008) is a visible photon counting detector designed as a 
guest facility on the SALT to provide very high time resolu- 
tion (<25 nanoseconds) and high signal to noise, full imaging 
photometry. UZ For was observed during a BVIT commis- 
sioning run on 25th January 2009 simultaneous with the 
HIPPO on the 1.9m telescope of the South African Astro- 
nomical Observatory. The data were extracted making use 
of the IDL data reduction software developed by the instru- 
ment team and binned into 0.5s bins. One eclipse of high 
time resolution and s/n was observed from which measure- 
ments were made. 



2.10 HIPPO January 2009, September 2010, October 
2010 and November 2010 

Unfiltered photo-polarimetric observations were made with 
the HIPPO (HI speed Photo-POlarimeter: Potter et al. 
2010) on four separate occasions and reduced as in Potter 
ct al. (2010). The single eclipse in January 2009 was ob- 
served simultaneously with the BVIT observations. Multi- 
ple eclipses were observed on the other occasions which were 
folded and binned from which measurements were made. 



3 RESULTS 
3.1 The eclipses 

A sample of our new and archival eclipse data is shown in 
Fig 1 phased on our new ephemeris (see below and table 2). 
The eclipse profiles are of varying quality and at multiple 
wavelengths. All of the eclipses can be understood in the 
framework of the standard polar model and from the gen- 
eral literature on UZ For (see section 1). UZ For undergoes 
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Figure 1. A sample of our new eclipse observations phased on our new ephemeris (section 4.2) Vertical grey bars indicate ingress and 
egress of the main accretion spot. Solid vertical bars indicate times of white dwarf ingress and egress assuming a duration of 40s. 



periods of either one (e.g. Bailey & Cropper 1991) or two- 
pole (e.g. Perryman et al. 2001) accretion states, which has 
been most clearly captured by our SALTICAM (2007) and 
BVIT (2009) observations respectively. During both accre- 
tion states the rapid fall and rise in flux at phases ~-0.031 



and ~0.031 correspond to the ingress and egress of the main 
accretion region. During the two-pole accretion state a sec- 
ond accretion region is additionally present, seen as the fall 
and rise in flux at phases ~-0.027 and ~0.027. The remain- 
ing gradual fall and rise during phases ~ —0.031 to —0.025 
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Table 1. Mid-eclipse times of the main accretion spot of UZ For. BJDyjjB is the Barycentric Julian Date in the baryccntric dynamical 
time system. Times have also been barycentrically corrected. 1 de Bruijne et al. (2002), 2 Perryman et al. (2001), 3 Imamura & Steiman- 
Camcron (1998), 4 Warren et al. (1995), 5 Ramsay (1994), 6 Bailey & Cropper (1991), 7 Allcn et al. (1989), 8 Ferrario (1989), 9 Bcucrmann 
et al. (1988), "'Osborne et al (1988). 
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and ~ 0.025 to 0.031 is attributed to the ingress and egress 
respectively of the white dwarf photosphere and takes about 
40s each. 



3.2 The O-C 

In table 1 we list all of our new mid-eclipse times as well as 
those we have measured from archival data or extracted from 
the literature. The orbital period calculation of Perryman et 
al. (2001) was used to calculate the cycle number for each 
eclipse. The period is sufficiently accurate to unambiguously 
assign cycle counts to the entire ~ 27 years of eclipses. The 
eclipse observed in our 2005 UCTPOL photometry was used 



to define the epoch (cycle 0). This period and epoch were 
next used as the starting point to perform a least-square 
quadratic fit with appropriate weighting set by the eclipse 
error measurements. The resulting fit gives a reduced \ 2 > 
95 with peak-to-peak residuals of ~ 60 — 80s. Note that 
we have not included the eclipses of Dai et al. (2010) in 
our analysis as their O-Cs are over 300s compared to our 
quadratic ephemeris. We believe either their measurement 
and /or time standard conversion to be in error. 

It is immediately apparent that there are significant de- 
partures from the quadratic ephemeris with a trend that ap- 
pears to be periodic (see top plot of fig 2 for residuals, albeit 
from a different quadratic fit: see below). 
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Table 2. Mid-cclipsc ephemerides of the main accretion spot of UZ For and corresponding planet model parameters. Ephemeris pa- 
rameters correspond to the representative solution (Fig 3) and are rounded off to the 1 sigma errors. The planet parameter errors were 
calculated using the range in parameter space of possible solutions and not the smaller 1 sigma errors of any one fit. Minimum planet 
masses are listed assuming coplancarity. M 3i j nc is the mass function. The combined mass of the primary and secondary stars are 
assumed to be O.84M0. 



Quadratic term: 



T = 2453405.30086(3) d 
P bi „ = 0.087865425(2) d 
A = -7(2)10- 14 



Planet 
Parameters: 



1st Elliptical term: 



v 3 = (E+T 3 )/ 3 

T 3 = 60383(416) (binary cycle) 

fz = 0.000098(3) (cycles/binary cycle) 

OT 3 = 0.85(5) 

K bin ,(3) = 0.00025(2) d 

e = 0.04(5) 



M 3}fnc = 2.9(1. 1)1O~ 7 M 
M 3tJup = 6.3(1.5) 
P 3 = 16(3) years 
a 3 = 5.9(1.4) AU 
oi, 2 = 0.042(1) AU 



2nd Elliptical term: 



v 4 = (E+T 4 )/ 4 

T 4 = 4833(215) (binary cycle) 
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K, 
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P A = 5.25(25) years 
a 4 = 2.8(5) AU 
01,2 = 0.025(1) AU 



c = 0.05(5) 



1983.7 1986.1 1988.8 1991.7 1993.9 1995.9 2000.0 2002.7 2005.2 2007.9 2010.9 

50 p . i ■ . . i . . . i . ■ ■ i ■ ■ . i . . . i — 




-8x10 4 -6x10 4 -4x10 4 -2x10 4 2x10 4 

cycle 



Figure 2. The O-C after successive subtraction of the three terms comprising our new eclipse ephemeris. Top: O-C after subtraction of 
the quadratic term with the first elliptical term overplotted (dashed curve). Middle: O-C after subtraction of the first elliptical term with 
the second elliptical term overplotted (dashed curve). Bottom: The final O-C residuals after subtraction of the second elliptical term. 
Diamonds are our new data or data that we have reduced from archives. Crosses are eclipse times from the literature and converted (by 
us) to BJDyjjs. 
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We next investigated the solutions resulting from mod- 
els consisting of a quadratic plus an elliptical fit to the eclipse 
times. A best reduced \ 2 °f 6.2 was achieved. An F-test 
shows that it is the better model (compared to the quadratic 
ephemeris) with a confidence of > 99.999% even though the 
elliptical term adds 5 more parameters to the model. 

However significant residuals still remain, (~10s) for 
some of the eclipse times (not shown). We next attempted a 
simultaneous quadratic plus two ellipticals fit to the eclipse 
times. The second elliptical term adds a further five param- 
eters to the model, giving 13 in total. Therefore, given the 
large number of parameters, a grid of starting parameters for 
minimisation was required in order to explore the resulting 
degeneracy in the solutions. Approximately 10 7 minimisa- 
tions were calculated. During minimisation all 13 parame- 
ters were free to vary. Predictably, the results have better 
reduced \ 2 but with degeneracy in many of the parameters. 
Formally the F-test confirms that adding a second ellipti- 
cal term is the better model with a high level of confidence 
(> 99.9999%) for the solutions with reduced \ 2 = 1-0. Other 
solutions with reduced \ 2 = 4-0 and 3.0 are also significantly 
better with a 98% and 99.9995% confidence respectively. We 
explore the degeneracy in the multi-dimensional x 2 space in 
section 4.2. 




1.5 1.6 1.7 1.8 1.9 2.0 

Frequency (1 CT 4 cycles/day) 



Figure 3. The reduced x 2 parameter space for the two elliptical 
periods. Black crosses, dark grey triangles and light grey dia- 
monds are the solutions with reduced x 2 < 1.0,1.0 < x 2 < 2.5 
and 1.0 < x 2 < 2.5 respectively. The diamonds have the ad- 
ditional constraint of both eccentricities < 0.1. The large black 
cross represents the location of the solution shown in Fig 2 with 
parameters listed in table 2. From left to right, the diagonal lines 
represent contours of constant period ratios of 3.1, 3.0 and 2.9. 
Typical one sigma errors are shown in the top left and bottom 
right corners. 



4 DISCUSSION 

4.1 The O-C 

Our results suggest that the deviations in the eclipse O-C 
are best described by the combination of a quadratic term 
plus two elliptical terms. This is highly suggestive of both 
secular and cyclic period variations. 

Period changes in binary systems are generally under- 
stood to be due to gravitational radiation, magnetic brak- 
ing, solar-type magnetic cycles in the secondary star (Ap- 
plegate's mechanism) and/or the presence of a third body 
in an orbit around the binary. 

Applegate's mechanism and/or the presence of a third 
body would be more consistent with cyclic variability. The 
latter would produce strictly periodic cycles while non- 
strictly periodic cycles would be expected from the former 
mechanism. Therefore, we next look at each of these mech- 
anisms in turn. 

4.2 Tertiary and quaternary components 

We now explore the degeneracy in the multi-dimensional 
X 2 space of the model fits containing one quadratic plus 
two elliptical terms in the context that the variations are 
due to the effect of third and fourth bodies in the system. 
Then the changes in the O-C arise because of the light-time 
effect caused by the gravitational influence of the additional 
bodies. 

As a first step, we plot the distribution of the period 
of the two elliptical orbits for those solutions which had re- 
duced x 2 < 2.5 (Fig 3). The starting grid, for minimisation, 
had period values between 2 and 50 years and Fig 3 shows 
that the minimised solutions have clustered in the period 
ranges ~ 13 — 19 and ~ 5 — 5.5 years. 

All of the solutions with reduced % 2 < 1-0 (x 2 < 28, 



shown as the black crosses in Fig 3) show periods centered 
on approximately 5.3 and 15 years, giving a period ratio 
of Rp ~ 2.8 . The corresponding predicted eccentricities of 
these solutions are >> 0.1 for both ellipses. We performed 
N-body simulations on a sample of these solutions. The Eu- 
ler method was used with sufficiently small time steps to 
ensure accurate calculations. We tested the accuracy of our 
code by first applying it to single elliptical orbits comparable 
to the innermost high eccentricity orbit. Our code preserves 
the eccentricity and the semi-major axis to better than 10%, 
and the periastron angle to < 0.1 radians over a time period 
corresponding to > 10 5 orbits of the outermost body. We 
then defined orbital solutions that have essentially the same 
two planet orbital elements to the starting conditions after 
> 10 5 orbits of the longer period as stable. As expected, we 
found that they all are very unstable orbits and are therefore 
unrealistic solutions. 

We next looked at the solutions with reduced x 2 m the 
range 1.0 < \ 2 < 2.5, which are shown as the dark grey 
triangles in Fig 3. They occupy a larger area of the plot 
which overlaps with the previous solutions. However most 
of these solutions still require the longer period elliptical 
to have a large eccentricity (> 0.1 and typically 0.4). N- 
body simulations within this parameter range also revealed 
unstable orbits. 

We therefore identified those solutions which had or- 
bital eccentricities < 0.1 for both planets and which had 
stable orbits according to our N-body simulations. These 
are represented as the light grey diamonds in Fig 3 which 
also overlap with the higher eccentric solutions presented 
above. However, the best of these solutions give X 2= 58 (28 
dof, giving x 2 =2.06), which is a significantly poorer fit com- 
pared to our best-fit highly eccentric orbital solutions and 
also gives a formally poor fit to our data. 

We note that our models assume the two planets to be 
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co-planar. There may yet be a set of realistic, more eccen- 
tric, solutions if the planets have inclined orbits with respect 
to each other. We have not investigated this additional pa- 
rameter space as our data set is not of sufficient quality nor 
quantity to warrant it. 

Given these caveats, in order to calculate the implied 
two planet parameters, we selected the best-fit low eccentric 
and stable solutions. Additionally we applied period errors 
that encompass the whole range of solutions in Fig 3. We 
note that the two planet parameters calculated below are not 
specific to this one best-fit solution but are representative 
of all the solutions with reduced \ 2 < 2.5 shown in Fig 3. 
The calculations are independent of the eccentricities. This 
particular solution is marked as the large black cross in Fig 3, 
the two elliptical parameters are listed in table 2 and over- 
plotted on the O-C in Fig 2. 

The amplitudes of the oscillations can be used to calcu- 
late the projected distances asin(i) from the center of mass 
of the binary to the center of mass of each of the triple sys- 
tems (0.042(f) and 0.025(1) AU for the long and short re- 
spectively). Setting the binary mass to be O.7M + O.14M 
gives the corresponding mass functions (/(m3,4) = 2.9 x 
fO~ 7 M ,5.3 x fO~ 7 M ). With the binary inclination at 
i = 80° the respective minimum masses for the third and 
fourth bodies (assuming they are in the plane of the bi- 
nary) are 0.006(f) M Q and 0.007(f) M Q and would there- 
fore qualify as extrasolar giant planets (6.3(1.5)Mj up and 
7.7(f.2)Mj„ p ) for orbital inclinations i 3 > 25°, u > 32° re- 
spectively. The quoted errors include the range in periods 
shown in Fig 3 and not the formal one sigma errors of one 
of the solutions. In addition, the quoted errors include the 
propagated uncertainties in the inclination and binary mass. 
These parameters are summarised in table 2. 

The equation for the times (T) of mid-eclipse of the 
main accretion spot are then given by: 

T(BJDtdb) = T + P bin E + AE 2 

(! - e (3)) 

+ ^n,(3)«n(«3-m 3 ) (1 + e(3)Cofl(t>3) 
+ ^n,(4)^4- CT 4) (1 + e(4)CoaM 

To, Pbin, A, E are the time of epoch, the binary orbital 
period (days), the quadratic parameter (related to the rate 
of period decrease by Pu n = 2A/P bin ) and the binary cycle 
number which comprise the quadratic term of the ephemeris. 
In the context that the two elliptical terms are due to third 
and fourth bodies in the system, then the parameters of 
the elliptical terms are: K b i n , ( 3j4 ) are the amplitudes of the 
eclipse time variations as a result of the light-travel-time ef- 
fect of the two bodies, W(3,4) are the true anomalies of the two 
bodies, which progresses through 2n over the orbital periods 
(P(3 j4 )) and are functions of E, the times of the periastron 
passages (T^^)) and the orbital frequencies (in cycles per 
binary orbital cycle) of the two bodies (/3,4). e<3,4) are the 
eccentricities and ^(3,4) are the longitudes of periastron pas- 
sage measured from the ascending node in the plane of the 
sky. Similar elliptical variations have been seen in NN Ser 
and DP Leo (Beuermann et al 2010a,b). 

Fig 2 shows the fit and the O-C residuals after successive 
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Figure 4. Solid curve shows the energy required to effect the 
period change observed in UZ For as a function of assumed shell 
mass, using Applegate's (1992) mechanism. The two horizontal 
lines represent the total radiant energy of the secondary (assum- 
ing 2880 < T e ff < 3020K) and hence the amount of energy 
available. 

subtraction of the three terms comprising our new eclipse 
ephemeris with parameters listed in table 2. The top and 
middle plots show that the two elliptical terms describe the 
time of eclipse variations very well. The lower plot shows the 
final O-C residuals after subtraction of the full ephemeris. 
Some residuals still remain which could be reduced further 
if larger eccentricities were permitted, particularly for the 
outer planet. Such orbital solutions may exist, especially 
given the possible indication that the planets could be locked 
in a 3:1 ratio. However better sampled observations are re- 
quired to further constrain the solutions. 

4.3 The secular variability 

The secular variability amounts to a decrease in the orbital 
period of P bi „ = -1.56(5) 10~ 12 s s _1 . Of the 1416 solutions 
that comprise the 'chosen' parameter space (light grey area 
in Fig 3) only one solution showed a Pbin consistent with 
0. It has a reduced \ 2 = 2-6. The rest predict a minimum 
in the period decrease rate of Pu n = —1.0 f0~ 12 s s~ 1 . For- 
mally an F-test shows that adding the quadratic parameter 
to the ephemeris is the better model with a 99.99% level of 
confidence (\ 2 = 109, 86 for 30, 29 degrees of freedom for 
the two models respectively). 

Similar levels of period decrease has been detected in 
other similar short period binaries e.g. DP Leo (Schwope 
et al. 2002), NN Ser (Brinkworth et al. 2006) and HU Aqr 
(Schwarz et al. 2009) for which gravitational radiation and 
magnetic braking have been shown to be either insufficient 
or problematic in the framework of the standard CV evo- 
lutionary model. Future observations are needed to show if 
these variations are indeed secular or periodic. 

4.4 Applegate mechanism 

Applegate (1992) proposed that solar-like magnetic cycles 
would drive shape changes in the secondary, thus redistribut- 
ing the angular momentum within the star, changing its 
oblateness. This then leads to a change in its quadrupole 
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moment and consequently a change in orbital period at the 
expense of some energy. This has been the preferred mecha- 
nism to explain such cyclic variations in CVs and other long 
period close binaries (Algol, RS CVn and W UMa stars). 

Following the prescription of Applegate (1992), the en- 
ergy required to generate a period change is: 

AJ 2 

AE = n dr AJ+— — 

21 eff 

The initial differential rotation f^r-AJ is set to zero in or- 
der to calculate the minimum energy required. The effective 
moment of inertia 7 e // = IsI*/(Is + I*) is calculated by 
assuming the star is separated into a shell Is and a core 7*. 
We experimented with a range of shell masses. AJ is the 
change in angular momentum and is given by 

-GMVayAP 
R \RJ 6tt 

We used the secondary star mass from Bailey & Cropper 
(1991) namely M = O.14M and the corresponding radius 
R = O.177i? following Patterson (1984). a = 5.5 x 10 8 m is 
the binary separation using a = 3.53 x 10 10 (Mi/M ) 1/3 (1 + 

l) 1/Sp olb( h ) (Warner 1995, equation 2.1b). AP can be ob- 
tained from equation (38) of Applegate (1992) relating the 
amplitude of orbital period modulation and the amplitude 
of the O-C oscillation: 

AP O-C 

— — — = Z7T— — 

where P and P mo d are the orbital and modulation period 
respectively (using Ps listed in table 2). The solid curve 
in Fig. 4 shows the minimum energy required to drive the 
maximum observed period change in UZ For, as a func- 
tion of assumed secondary shell mass. The two horizontal 
lines represent the total radiant energy of the secondary 
over the modulation period, assuming 2880 < 
T e ff < 3020K, which appears to be more than sufficient to 
drive the Applegate mechanism. The situation is not so clear 
cut if one instead integrates over shells and allows for the 
quadrupole moment of the core (using the calculations of 
Brinkworth et al. 2006). This raises the minimum energy by 
about an order of magnitude, which makes it comparable, 
at minimum, to the energy of the star. 

We should add that Lanza et al. (1998) propose a pre- 
scription that is more energy efficient than the Applegate 
mechanism, perhaps by a factor of two. Therefore, with fur- 
ther refinements, magnetic fields may yet be shown to able 
to drive the period changes seen here. 



4.5 Spot motion 

The eclipse times are derived from the observed ingress and 
egress times of the accretion spot and not the center of the 
white dwarf itself. Therefore the observed O-Cs could be as 
a result of motion of the spot on the white dwarf. In addi- 
tion, the egress and ingress of the white dwarf photosphere 
has been observed to take about 40 seconds (seen unam- 
biguously in the low state observations of Bailey & Cropper 
(1991) and confirmed in our low state SALTICAM 2007 ob- 
servations) and therefore could accommodate a ~ 40 — 60s 
of spot motion. 



To assess this possibility we investigated the actual mor- 
phologies of the eclipse profiles in order to measure the rel- 
ative phases of the white dwarf photosphere ingresses and 
egresses to that of the main accretion spot. Accordingly, we 
display a sample of our eclipse observations in Fig. 1 phase 
folded and binned on our new ephemeris. The upper two 
eclipse profiles are the HIPPO 2010 and the BVIT 2009 ob- 
servations respectively. These correspond to the two newest 
data points in Fig. 2 which show a O-C shift of ~ 30s with 
respect to each other on the quadratic subtracted O-C plot 
(upper plot Fig. 2). Therefore if the spot did indeed change 
locations, on the surface of the white dwarf, between these 
two observations then we should observe a corresponding rel- 
ative phase shift between spot ingress/egress to that of the 
white dwarf photosphere ingress/egress. However one can 
see that the beginning of the white dwarf photosphere egress 
(solid vertical line at phase ~ 0.026) is consistently ~ 40s 
(0.005 phase) ahead of the spot egress (vertical grey bars) 
between these two observations. Thus, one would have ex- 
pected the relative time difference between the white dwarf 
photosphere and spot to be shorter by ~ 30s between the 
two observations and not approximately equal as observed. 
Therefore we conclude that the spot has not moved on the 
surface of the white dwarf during these observations, at least 
within our measurable errors of about l-5s. The same rela- 
tive phase difference (spot, white dwarf photosphere) is also 
apparent in the other eclipse profiles in which the spot and 
white dwarf photosphere are resolved (see the next five plots 
in Fig. 1 corresponding to the SALTICAM 2007, UCTPOL 
2005, SWIFT 2002, UCTPOL 2002 and SCAM 1999 obser- 
vations). Furthermore, the same unchanging relative time 
differences are seen in the ingresses, although the white 
dwarf photosphere ingress emission may be complicated by 
an additional contribution from the accretion stream. The 
stream is not visible during the white dwarf photosphere 
egress, which can be understood from simple eclipse geo- 
metrical arguments, and confirmed through HST UV spec- 
troscopy (Stockman & Schmidt 1996). 

We note, however, that the longitude of the accretion 
spot should be expected to change during different accre- 
tion states. For example, Schwope et al. (2001) calculated a 
change in spot longitude of ~ 10° between high and inter- 
mediate accretion states in the eclipsing polar HU Aqr. This 
would translate to about a shift of ~ 2 — 3s in the O-C values 
(Schwarz et al. 2009) which therefore cannot account for the 
observed O-C values. This implies that if there was a similar 
spot motion in UZ For during different accretion states, it 
cannot account for the large shift seen in the measured O-C 
values. 

Additionally there was not any measurable movement of 
the spot in latitude during our observations. This is apparent 
from table 1 where the eclipse width measurements agree 
within errors: a change in spot latitude would have resulted 
in a corresponding change in eclipse length. 



5 SUMMARY AND CONCLUSION 

We have detected departures in the eclipse times of UZ For 
from a simple quadratic ephemeris of up to ~60s. The de- 
partures are suggestive of two periodicities of ~16 and ~5.25 
years. The two favoured mechanisms to drive the periodic- 
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ities are either two giant extrasolar planets as companions 
to the binary or a magnetic cycle mechanism (e.g. Apple- 
gate's mechanism) of the secondary star. However, Apple- 
gate's mechanism would require the entire radiant energy 
output of the secondary and therefore would seem to be the 
least likely of the two, barring any further refinements in 
the effect of magnetic fields (e.g. Lanza 1998). A two planet 
model is also problematic given the quality of the data in 
that a high eccentric orbit, for at least one of the planets, 
seems to be required to fully capture all of the eclipse times. 

If it can be confirmed that the residuals are due to a 
third and a fourth body, then the planets either formed in 
a pre-common envelope circumbinary protoplanetary disc 
(first generation) or in a disc that resulted from the com- 
mon envelope (CE) phase (second generation: Perets 2010). 
The separation of the progenitor binary is of the order of a 
few AU, comparable to that of the planets, which implies 
that only second generation planets could have formed at 
the orbits suggested here. However, Beuermann et al (2010a) 
suggested for the planets around NN Ser, a slowly expanding 
CE could provide the dynamical force to drag inwards plan- 
ets formed further out, which would have otherwise been 
lost to the system due to the decrease in mass of the central 
binary (Alexander et al. 1976). In either case, we note that 
the semi-major axis of even the shortest period object poses 
no problem for orbit stability (Holman & Wiegert 1999). 

It is intriguing that Qian et al (2011) propose a very 
similar two elliptical model fit for the polar HU Aqr, also 
using eclipse timing results. In particular they also find that 
the larger ellipse requires a high eccentricity (0.51) to cor- 
rectly capture all of the data. Therefore, their two planet 
model for HU Aqr seems to have a similar problem with 
orbit instabilities that we have found for UZ For. 

As yet there is insufficient data on UZ For to iden- 
tify conclusively the mechanism responsible for the periodic 
changes in its eclipse times, and indeed more than one mech- 
anism could be present. Further good signal/noise, high time 
resolved observations of UZ For and other similarly eclipsing 
systems are encouraged. 
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